Influence of the pseudogap on the superconductivity-induced phonon renormahzation 

in high-Tc superconductors 
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We investigate the influence of a d-density wave (DDW) gap on the superconductivity-induced 
renormalization of phonon frequency and linewidth. The results are discussed with respect to 
Raman and inelastic neutron scattering experiments. It turns out that the DDW gap can enhance 
the range of frequencies for g = phonon softening depending on the underlying band structure. 
Moreover we show that an anisotropic 'd-wave' pseudogap can also contribute to the q-dependent 
linewidth broadening of the 340cm~^ phonon in YBa2Cu307. 



The origin of the pseudogap is presently a strongly 
debated issue in the field of high-Tc superconductivity 
and has been detected by numerous techniques in all the 
cuprates 0]. One of the ideas which received consider- 
able attention involves the existence of preformed pairs 
This idea finds support in continuous evolution of 
angle-resolved photoemission spectra (ARPES) from the 
normal to the superconducting (SC) state However, 
another appealing scenario suggests that the pseudogap 
could (not only) arise from pairing in the particle-particle 
channel, but also from different scattering mechanisms 
(like CDW fluctuations) which become singular near op- 
timal doping. Indeed there is increasing experimental ev- 
idence that the peculiar properties of the cuprates, 
both in the normal and the SC phase are related to the 
occurrence of a Quantum Critical Point (QCP) located 
near (slightly above) the optimal charge carrier concen- 
tration. Besides singular scattering at the QCP a further 
direct consequence of this scenario is the prediction of an 
ordered state in the underdoped regime. Various types of 
orderings have been proposed each of them breaking dif- 
ferent kinds of symmetries. For example, Castellani et al. 
Q have proposed incommensurate charge-density waves 
(ICDW) which break translational invariance. Varma 
[0 considers circulating currents which (besides time- 
reversal symmetry) break a rotational invariance of the 
copper oxygen lattice but preserve translational invari- 
ance. In one of the latest proposals Chakravarty et al. 

have investigated d-density wave (DDW) order. This 
state is made of staggered currents which break parity 
and time-reversal symmetry as well as translational in- 
variance by one lattice constant and rotation by 7r/2. All 
of these proposals have in common that the correspond- 
ing order parameter induces a gap in the single-particle 
spectrum which has the same anisotropy as observed in 
ARPES experiments. For the circulating current state 
this has been shown in Rcf. |^ and for the ICDW phase 
m Ref. 0. In case of the DDW state the pseudogap 
anisotropy is explicitly put in the quasiparticle hamilto- 
nian. 

In this paper we consider the influence of an anisotropic 
single-particle gap on the frequency shift and linewidth 



of phonons in the SC state. Strong phonon self-energy ef- 
fects in high-Tc cuprates have been analyzed first by Zey- 
her and Zwicknagel within Eliashberg theory. For an 
isotropic s-wave superconductor they obtained softening 
for phonons with frequencies uj < 2A and hardening for 
phonons with uj > 2A. In addition the linewidth of the 
phonon mode is only affected for uj > 2 A where it is addi- 
tionally broadened due to the opening of the SC (s-wave) 
gap. The Zey her- Zwicknagel approach has been general- 
ized to anisotropic gap functions and dispersions relevant 
for the cuprates in Refs. |T2|jr^ ] . The main difference with 
respect to the s-wave model is the appearance of an ad- 
ditional line broadening for mode frequencies to < 2Aq 
since in a d-wave superconductor the density of states 
(DOS) remains finite for energies within the maximum 
SC gap. 

Our investigations are based on the DDW state since 
from the formal point of view this scenario offers a simple 
way of introducing an anisotropic gap into the problem. 
However, we want to stress that most of the results pre- 
sented below should also (at least qualitatively) hold for 
other QCP proposals as long as the corresponding pseu- 
dogap leads to an analogous quasiparticle dispersion. In 
this context it is interesting to note that a DDW-type 
gap has also been considered for the ICDW scenario . 
Starting point is the following hamiltonian 
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where x(k) = Xo[cos{kx) — cos(fcj,)]/2 and A(k) — 
Ao[cos(fca;) — cos(fcj^)]/2 denote the DDW and SC gap re- 
spectively and Q = (tt, tt). The dispersion will be spec- 



ified later. In the basis = (c^i, c_kj,, c^+qp c. 
the hamiltonian eq. (|^) is given by 
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where Sk = (sfe - £fe+Q)/2 and fl^ = ^j, - {sk 
Diagonalization yields the four eigenvalues 
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FIG. 1. Density of states in the superconducting state (Ao — 
23meV) for three values of pseudogap: xo/^o = - solid 
curve,xo/Ao = 1 - dashed curve, xo/Aq — 2 - dotted curve, 
a) Nearest-neighbor tight-binding model at half filling; b) En- 
ergy dispersion from Ref. M]. 
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A| , where Ek 



the reduced zone of the DDW state g . 

The superconductivity induced phonon renornial- 
ization is obtained from the imaginary-time or- 
dered charge-charge correlation function I{{q,iuj) — 
{TrPq(iuj)p-q{—iw)) via = (7^n(q, w-|-i(5) and the cor- 
relator is analytically continued to real frequencies at the 
end of the calculation. Throughout the paper we restrict 
ourselves to a momentum independent electron-phonon 
coupling constant g in order to focus on the essential 
contribution from the pseudogap anisotropy. 

We start our analysis by calculating the renormaliza- 
tion of g = phonons. In this case the charge suscepti- 
bility reads as 
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where ns{k) = [e^^^'^'^'i + 1)^"'^. Note that within our 
approach the q=0 phonon self-energy is only finite in the 
SC state for arbitrary value of the pseudogap. In order to 
investigate the influence of the pseudogap on the phonon 
renormalization above Tc one should include vertex cor- 
rections to n(g — 0, iuj) which however is beyond the 
scope of the present paper. 

Obviously the underlying dispersion plays a significant 
role in determining the properties of T\{q = OjCj). Fig. 
1 shows the density of states (DOS) for two model band 



dispersions which will be considered below. The DOS 
is evaluated in the SC state (Aq = 23meV for all fol- 
lowing results) and for three values of the pseudogap. 
Fig. la corresponds to = —2t[cos{kx) + cos{ky)] with 
t = O.beV i.e. complete nesting for the DDW scattering 
vector at half-filling. Due to particle-hole symmetry the 
chemical potential fj, is not altered in this case upon in- 
creasing the pseudogap. Thus for ^ = xi^) ^tnd A{k) 
work cooperatively and the DOS shows a single 'lead- 
ing edge' gap at energy 2-\/A^ + Xo- ^ig. lb shows the 
DOS for a tight-binding dispersion recently proposed 
by Norman in order to analyze neutron-scattering 
and ARPES data in the SC state of underdoped bilayer 
cuprates. In this case the gap edges of SC gap and pseu- 
dogap are separated since the former opens at the chem- 
ical potential whereas the pseudogap is related to the 
energy where the dispersion is most susceptible to DDW 
scattering. Moreover the DDW state induces a Fermi 
surface which consists of four pockets centered around 
(±7r/2, ±7r/2). With increasing DDW gap the pockets 
shrink towards the nodes thus decreasing the effective 
SC gap in the DOS (cf. Fig. lb). 

For both dispersions Fig. 2 depicts the real and imag- 
inary part of n((7 — 0,uj) as a function of cij/(2Ao) and 
for zero temperature. In case of zero pseudogap the fre- 
quency shift and linewidth broadening correspond to the 
results of Refs. [|l2|,^ i.e. crossover from hardening to 
softening at frequencies uj « 2Ao. In addition the d-wave 
gap results in a linewidth broading also for oj < 2A due 
to the finite density of states in the SC gap. However, 
consider the case of finite DDW gap and the underly- 
ing dispersion (left panel). It turns out that the sign 



change in Slo now occurs at frequencies 2-\/Aq 
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FIG. 2. Real and imaginary part of n(g = 0, uj) corresponding 
to the DOS shown in Fig. la (left panel) and Fig. lb (right 
panel). Labeling of the curves is the same than in Fig.l. 
The dashed lines in the upper right plot indicate the phonon 
frequencies discussed in the text. 
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following the behavior in the DOS. Moreover with 
increasing pseudogap the phonon linewidth broadening 
rapidly decreases and correspondingly also the maximum 
phonon frequency shift. In case of the dispersion e'j. (right 
panel) the situation is more complex and can be analyzed 
best by first considering the linewidth broadening which 
follows closely the behavior of the DOS. Upon increas- 
ing the DDW gap Imll{q — 0, w) develops a double peak 
structure where the low energy excitation is due to 
SC and the other one is due to the DDW 'leading edge' 
in the DOS. This structure is reflected in Rell{q = Q,uj) 
by the two step like features below and above 2A(0). 
It turns out that at least up to x/^(0) = 1 the DDW 
gap enlarges the frequency range where phonon soften- 
ing occurs. In addition the 'hardening peak' becomes 
strongly suppressed by a finite value of the pseudogap 
most notably exemplified in case of x/A(0) = 2 where 
Rell{q — 0,uj) is practically zero above oj ^ 2A(0). 

We now turn to the analysis of the temperature de- 
pendence of the phonon renormalization which we eval- 
uate for the dispersion e^. Fig. 3 shows the real and 
imaginary part of n(g = 0, lu) for the two frequen- 
cies uja,b = 2A(0) ± A(0)/2 (indicated by dashed lines 
in Fig. 2b). For simplicity, we adopt the BCS tem- 
perature dependence for both SC and DDW order pa- 
rameters, i.e. Ao(T) = ^1 - (r/Tc)2,Tc = 80K and 
Xo{T) = y/1 - {T/T*Y,T* = 150/^: (note that the tem- 
perature dependence of xo {T) has no significant influence 
on the result). For zero DDW gap the a-phonon initially 
hardens close to Tc where the 'effective' SC gap is well 
below uja in agreement with Ref . JTst . This hardening be- 
comes completely suppressed upon increasing the DDW 
gap, however, there is still significant softening for T ^ 0. 
Correspondingly the peak in Imli{q — 0,a;) is also re- 
duced by the opening of the pseudogap. As already dis- 
cussed above the hardening of the b-phonon is much more 
affected by the DDW gap. In fact the frequency shift is 
almost completely suppressed for x/A(0) = 2 where also 
the strongly reduced DOS in the pseudogap reduces the 
phonon scattering and thus the corresponding linewidth. 

Our results can be compared with doping dependent 
Raman scattering experiments where upon underdoping 
the increasing pseudogap should have a significant in- 
fluence on the phonon renormalization. The most pro- 
nounced features are found for the 340 cm~^ phonon in 
YBa2Cu307_2: which undergoes a strong softening and 
hardening of the 440cm~^ and 500cm~^ modes below the 
SC transition temperature. The particular role of dop- 
ing for the Raman modes below Tc was investigated by 
Altendorf et al. [0 and more recently by Limonov et 
al. [Q. It turns out that the superconductivity induced 
changes in linewidth and frequency for the 340cm~^ Big 
mode depend very sensitively on the charge carrier con- 
centration near optimal doping. With increasing under- 
doping the phonon renormalization becomes rapidly 
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FIG. 3. Temperature dependence of frequency shift and 
linewidth broadening for two frequences LOafi = 2A(0) ± 
A(0)/2 (left and right panel). Labeling of the curves is the 
same than in Fig.l. 

suppressed in contrast to the overdoped regime where 
the data of Ref. [|l8| do not show a pronounced difference 
of frequency shift to the optimally doped case. Only the 
linewidth broadening gradually decreases from the over- 
doped to the optimally doped sample and vanishes for the 
underdoped system. Analagous behavior is observed for 
the 430cm~^ and 500cm~^ modes which show hardening 
for the optimal and overdoped samples while frequency 
shifts are suppressed in the underdoped regime. Thus 
these features can be (at least qualitatively) accounted 
for by evaluating the phonon renormalization with an 
anisotropic pseudogap as described above. 

We now proceed in calculating the influence of the 
DDW gap on the q-dependent phonon renormalization 
in the SC state (based on dispersion e\). For simplic- 
ity we restrict ourselves to zero temperature where the 
charge-charge correlation function is given by 



n(g, icj) = 



1 

2N 



E 



Es{k + q) + Et{k) 



{iujY ~ {E,{k + q) + Et{k)y 



The coherence factors read as 
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Since we are interested in the change of phonon fre- 
quency and linewidth upon entering the SC state we 
consider in the following the difference n'^>*'(q, icj) — 
Il^^^{q,iu!). Fig. 4 depicts the corresponding real and 
imaginary part for uia = 2A(0) — A/2 along the [1,0] — 
and [1, 1] direction for different values of x(0). Note that 
the peaklike structures at small q are due to the two- 
dimensionahty of the Fermi surface (see e.g. [pj|j20|) and 
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FIG. 4. SC induced change in the q-dependent frequency shift 
and linewidth broadening for cua ~ 2A(0) — A(0)/2. Labeling 
of the curves is the same than in Fig.l 



deduce that the measured linewidth broadening along 
the (1, 0) direction is at least partially due to the pres- 
ence of the pseudogap. This feature together with the 
reduced frequency renormalization of phononic Raman 
modes upon underdoping supports the point of view, that 
the pseudogap also persists below Tc and may not (only) 
originate from pair fluctuations. 

Obviously for a quantitative comparison a more de- 
tailed analysis for the q-dependent electron-phonon cou- 
pling is necessary which for the 340cm~^ mode has been 
done in Ref. ||2^. In addition, we have restricted our- 
selves to the purely two-dimensional BCS case. How- 
ever, at least the qualitative features of our investigations 
should be robust with respect to these approximations 
and may provide a basis for a more quantitative analy- 
sis of phonon renormalizations under the presence of an 
anisotropic pseudogap. 

The authors would like to thank A. Bill for helpful 
discussions and a critical reading of the manuscript. 



they are only observable when the condition Qz = is 
fulfilled. However, since the latter is experimentally hard 
to realize we disregard these features in our following 
analysis. 

Consider first the (0, 0) — > (tt, 0) (F M) scan where 
the pseudogap induced change in ReIV{q, oj) is most sig- 
nificant for wave vectors up to g ^ (7r/2,0). In fact, for 
zero DDW gap the system is rather susceptible for per- 
turbations within this wave vector range since the corre- 
sponding scattering acts in the high-density (tt, 0) region 
of the underlying dispersion. However, it is exactly this 
part of the Brillouin zone where the DDW state induces 
a gap which in turn leads to the decrease in the charge 
susceptiblity for < jl. As a consequence the most 
pronounced change in the SC induced phonon frequency 
shift now occurs at around q w (7r/2, 0) which reflects it- 
self in an increase of the linewidth around the same wave 
vector. For the (0, 0) (tt, tt) (F X) scan the same 
arguments than above can be applied in order to under- 
stand the suppression of small q frequency shifts upon 
increasing the DDW gap. Moreover, along the diago- 
nals also the large wave vector phonon renormalization 
is affected. Thus the opening of the pseudogap strongly 
reduces the SC induced frequency change over large parts 
of the F — > A" direction. 

Therefore, our results can account for some of the 
anomalies of the 340cm~^ phonon observed by inelas- 
tic neutron scattering [El 22|. Most interestingly it was 



found that at low temperature the SC induced change 
in phonon frequency is largest along the (1,0) direction 
for wave vectors up to g ~ (7r/2,0). Then the softening 
progressively vanishes towards the zone boundary and 
as a consequence the linewidth aquires a maximum at 
q « (7r/2,0). On the other hand much smaller soften- 
ing and no linewidth broadening has been detected along 
the diagonal direction. Thus from our analysis one may 
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